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Abstract 

In the reversible Schrodinger-Newton equation a complex Newton coupling 
Gexp(— za) is proposed in place of G. The equation becomes irreversible 
and all initial one-body states are expected to converge to solitonic station- 
ary states. This feature is verified numerically. For two-body solutions we 
point out that an effective Newtonian interaction is induced by the imaginary 
mean-fields as if they were real. The effective strength of such induced gravity 
depends on the local wave functions of the participating distant bodies. 
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1. Introduction 

The Schrodinger-Newton equation (SNE) was proposed in the context 
of quantum foundations [H, [ij. Assuming a Newtonian self-interaction, it 
realizes the concept that gravity is responsible for the observed spatial local- 
ization of macroscopic objects. It leads to plausible scales of localization. The 
SNE has been studied in numerous works @, @, 0, 0, @, 10, 11, 12, 13, 14, 15 



Its basic feature is that the Newtonian mean-field self-interaction modifies 
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the free Schrodinger equation in such a way that it acquires locahzed ground 
states (sohtons). These solutions are then considered as natural 'pointer' 
states for the macroscopic objects. On the other hand, these stationary 
pointer states are expected to emerge through the evolution of the wave 
function. Unfortunately, this irreversible mechanism can not be realized by 
the reversible SNE. To make the solutions converge toward pointer states, 
one gives the Newton coupling G a negative imaginary part. It turns the re- 
versible SNE into an irreversible frictional SNE (frSNE). The specific frSNE 
with G — > —iG has been an emergent structure in the author's decoherence 
model 0, 01, and its relationship to the reversible SNE has been analyzed 



in Refs. [15| and [16|. We recall the one-body frSNE in Sec. [21 we solve 
it numerically and determine the unique stationary wave packet in Sec. O 
We discuss the two-body frSNE in the special case when the bodies are far 
form each other and find that the imaginary coupling induces gravitational 
attraction, as if it were real coupling (Sec. H]). The strength of this induced 
gravity, however, turns out to depend on the quantum state of the partici- 
pating bodies. The general case of complex coupling Ge"*" with < a < vr 
is briefly discussed in Sec. O 

2. One body equation 

The SNE is a Schrodinger equation with a Newtonian mean-field potential 
where - in contrast to textbook many-body equations - we retain the self- 
interaction terms. Therefore the mean-field potential is already present in 
the one-body SNE: 

where ip = ip{r) is the c.o.m. wave function of the object of mass M, G is 
the Newton constant, and 

= V^r) = -GM' J y^^dr' (2) 

is the Newtonian mean-field potential of a point-like object. The stationary 
ground states are solitons of spread (Ar)o ~ {h^/GM"^) but they do not 
attract other solutions. We create a basin of attraction (convergence) if we 
give G an imaginary part. Robust convergence is achieved by the replacement 

^-^A,-i(V,-(nW. (3) 
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The imaginary mean-field corrupts the normahzation of we restore it by 
the constant counter term: 



{V,) = -GM- j '^^ifj^Ji^^'^drdr- . (4) 

The frSNE preserves both the momentum and the position expectation values 
(p) = p and (r) = r. If, for simphcity, we start a wave function with 
p = and r = then the solution will converge to the rotational invariant 
stationary state 

V^o(r)e-(^/^)^°* , (5) 



where 



Eo = ^ I mr)\U7rr'dr . (6) 



Surprisingly, the 'energy' Eq of the stationary state contains the kinetic 
energy only, without the contribution of the Newton self-interaction. This 
is a consequence of the choice G —iG which makes the mean-field pure 
imaginary, there is no real dynamical potential left in the frSNE. 

Little is known about the details of the stationary solution ([5]). Dimen- 
sional analysis of the standard spread yields the same order of magnitude 
(Ar)o ~ (h^/GM^) as in case of the reversible SNE Analytic results ex- 
ist for large extended spherical objects 0, [lij; for point-like bodies we have 
to use numeric simulations. Once we know the stationary solution i/jq in the 
c.o.m. frame then, thanks to the Galilean invariance of the frSNE, we know 
all stationary solutions: 

V^odr-r.De^^/'^^^P^-^*) , (7) 

where E = Eq + (pV2M) and dvt/dt = p/M while p = const. If the frSNE 
is perturbed by a weak smoothly varying field, real or imaginary, we can still 
retain the above form if we make the c.o.m. momentum p time-dependent. 
This will be the case later when we apply the frSNE to two bodies far from 
each other. The following position-momentum correlation matrix will play a 
definitive role in their effective interaction: 

1 —if 1_ _ 

R = ((P ° r) — (p) o (r)) = — / '?/'*(r)(V o r + r o S/)ip{r)dr — — p o r . 



(8) 
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Due to the rotational invariance of ipo, the stationary matrix Rq becomes 
proportional to the unit matrix I, i.e.: Rq = -Rq'- The correlation scalar Rq 
will be determined numerically. 

3. Numerical solution 

We restricted our numeric simulations for rotational invariant states ip{r) = 
ip{r). The mean-field takes the following form: 

VJr) = -AnGM' I '^^^f^' r'^dr'. (9) 
J max(r, r') 



The frSNE reduces to: 

We have simulated the solution of this equation with various initial wave 
functions (Gaussian, smoothened rectangle, superposition of two Gaussians). 
They all converged to the unique localized stationary solution We il- 
lustrate the convergence process by monitoring the time-dependence of the 
standard spread Ar of the wave function: 

oo 

Att f 

(Ar)2 = — |V^(r)|VMr . (11) 



A competition between the spreading kinetic term and the contractive po- 
tential term is witnessed by transient oscillations until the stationary value 
(Ar)o = 5.5501 (^V<^^^) is reached (Fig. [I]). The Fig. [2] shows the station- 
ary distribution |'?/^o('^)P- It is bell-shaped but it is not Gaussian. Also the 
complex phase of the numeric solution ^/'o(^) is shown in the same figure. 
We numerically observed the constant shift of the phase at speed —Eq/H, 
with Eq = 0.0356(G^M^//i^). This value has independently been confirmed 
by numeric integration of the expression ([6]) using the numerically calculated 
ipQ. We also determine the stationary position-momentum correlation matrix 
([8j) in the numerically obtained stationary state ipQ. The correlation scalar 
Rq can be written as: 

oo 

4-Tr f 

i?o = ^y r'mr)?^,{r)<lT , (12) 
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Figure 1: Relaxation of the standard spread Ar of the wave function ip{r,t) toward the 
stationary value (Ar)o = 5.5501 in function of time. A real Gaussian of spread 1 was the 
initial wave function [h = G = M = 1). 




Figure 2: Spatial density \ipo{r)\'^ (full) and phase Xoir) (dashed), resp., of the stationary 
wave function tpo{r) in function of the radial coordinate r. Standard position spread is 
(Ar)o = 5.5501, momentum spread is (Ap)o = 0.2668, position- momentum correlation is 
Ro = 0.6753 {h = G = M = 1). 
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where Xo = arg('?/'o)- Numeric integration obtained Rq = 0.6753. 



4. Two body equation 

For two identical objects the wave function reads \l/(ri, and the frSNE 
has the following structure jl, 111] : 

f = ^A,*-i(vr+v^-'-(<"+vr))* 

+^A,* - + " - (v;;-' + v'-i))* , (13) 

where the Laplacians Ai,A2 refer to ri,r2, respectively. The mean-field 
self-interactions read: 



l^(r'i,r^: 
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Vi''\r^) = -GM' / '^^^^'^^ dr^dr^ , (14) 
Vr\r2) = -GM' I 'T^^^'^V,' dr^dr; , (15) 



the mean-field cross-interactions read: 



|ri - 








|r2 - 




l^(r'i, 




|ri - 




l^(r'i, 




|r2 - 


ril 



Vj' ^(ri) = -GM' I j/,' dr;dr- , (16) 



K{"^(r2) = -GM^ / ^Aiililf dr;dr'2 . (17) 



We do not intend to discuss the generic two-body solutions. We are 
interested in the special case when the two bodies are far from each other 
at locations fi = (ri) and r2 = (r2). If |ri — r2| is much larger than the 
one-body stationary extension (Ar)o then we expect the following scenario. 
Considering first the limit |ri — r2| ^ oo, the cross-interactions and 
V} can be omitted, the two-body frSNE (|T3l) splits into two separate one- 



body equations ([3]): both bodies reach their localized stationary states in 
their respective c.o.m. frames. Accordingly, the two-body solution is simply 
the product 

vi/(ri,r2,t)=4'^(ri,t)4')(r2,t) (18) 
of the one-body solutions ([7]): 

4'\rut) = ^odn - ritDe^^/'^)^^^''-^^*) , (19) 
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where Ei = Eq + (p^/2M) and dru/dt = p^/M, and a similar definition 



(2) 

holds for tpQ , too. As we said, the effect of the cross-interaction terms 



(12) (21) 

is perturbative and leads to the slight acceleration of the c.o.m. 
momenta Pi,P2- Let us calculate the acceleration of p^^. First, we expand 

(12) 

the cross-interaction term . 

C^(ri) = C^(ri)-F(ri-fi), (20) 

where: _ 

I = CM' ~_^'3 . (21) 
|ri - r2| 

If we substitute the ansatz f|T8|) into the two-body frSNE f|T3l) . the following 
separate structure can be obtained for ipl^'^: 

This is an effective frSNE, valid for (Ar)o <^ |ri — f2|, to calculate the c.o.m. 
acceleration dp^/dt of the stationary wave packet of the body. The vector 
F is the Newton force at location ri = (ri) caused by the other body at the 
remote location f2 = (r2), and — F is the force on the second body caused 
by the first one. A similar equation could be derived for dpg/dt, with the 
opposite force — F. However, the vector F does not play the role of a real 
force in the Eq. (122|) . If it did, it should have come with the additional 
imaginary factor i. Nevertheless, we are going to prove that even these 
imaginary Newton forces can mimic the true Newton forces. Indeed, they do 
accelerate the bodies in a Newtonian way apart from a numeric factor Rq. 
The acceleration dp^/dt has the standard form: 



— = -ihJ [tp'o '{n)] Vi dri -thj 



d^i'\ri) 
dt 



Vi4'Vi)dri 



(23) 

Let us substitute the effective equation fl22|) and recall the matrix Rq for the 
position- momentum correlation ([8]) in the stationary state tpQ. After trivial 
steps we get dp^/dt = 2RoF which, due to the rotational symmetry, amounts 
to the ultimate form: 

^ = 2i?oF ^ 1.3506F , (24) 
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Figure 3: The ratio of the effective coupling Ga to Newton's G in function of the phase a 
of the complex Newton coupling Ge"*". 



and, of course, similar steps would give the opposite acceleration dpa/dt = 
—2RqF k, — 1.3506F for the other body. The value of Rq was obtained 
numerically in Sec. [3l 



5. Complex coupling 

We outline some features of the general frSNE - first advocated in the 
present Letter - where the coupling G in the SNE is replaced by Ge"*". 
Obviously, a = means the reversible SNE ([T]), a = 7r/2 is the case G —iG 
known earlier 0, 16 1 and studied through Sees. [2}|1] while a = vr would mean 



the reversible SNE with repulsive Newton gravity which we exclude from 
considerations. Inside the region < a < vr we have generalized frSNEs with 
features resembling the investigated special case a = Tr/2. All these frSNEs 
have localized one-body asymptotic states of unique shape with parameters of 
the same order of magnitude. Convergence was numerically observed already 
at small imaginary couplings. (We conjecture that the ground state solution 
of the SNE is exactly indetical with the stationary state ipQ of the frSNE 
with coupling G — ie, where e +0.) For a = tt/2 one has a fairly robust 
convergence which becomes finally lost towards the edges a = 0, vr. For a 
general frSNE, the strength of Newton attraction has two contributions: the 
true dynamics which is proportional to ReG and the induced acceleration 
which is proportional to 2RolmG. The resulting gravity is characterized by 
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an effective coupling: 



Ga = (cos a + 2Rq sin a)G . (25) 

On Fig. [31 we have plotted the factor cos Q; + 2i?o sin a in function of the phase 
a. We see that the effective Newton forces are bigger than the standard 
ones by a small numeric factor about, apparently, less than 2. We could 
renormalize the bare coupling G by the inverse of this factor in order to get 
the right effective coupling. This discussion is beyond our scope now, yet 
we mention a very similar issue with an alternative - and perhaps related - 



concept of emergent (induced) gravity 17 



6. Summary 

We have studied the simplest frictional Schrodinger-Newton equation 
characterized by the pure imaginary coupling —iG. We solved the one-body 
frSNE numerically and found robust convergence to a unique localized wave 
packet in the c.o.m. frame. In case of the two-body frSNE we discussed a 
heuristic solution in the special case when the two bodies are far from each 
other. They quickly form stationary wave packets in their respective c.o.m. 
frames. An effective attraction emerges between them which is 2Rq times 
the Newtonian attraction where Rq is the position-momentum correlation of 
the one-body stationary wave packet. Therefore our 'induced gravity' de- 
pends on the details of the wave function of the participating bodies. Before 
they reach their stationary states towards the limit 2R 2Rq ^ 1.3506, the 
emergent gravity may be very different from Newton's. As to the station- 
ary regime itself, an immediate question arises: can we tune the factor Rq 
to 1/2? Yes, we can. The frSNE of extended spherical bodies yields just 
Rq = 1/2 0, 17 1 . Therefore it may well be that the concept of the frSNE 



with pure imaginary coupling —iG contains a bit of real physics: for large 
extended objects the frSNE would show stable convergence to one-body lo- 
calized states and induces correct Newton forces between the distant objects. 
Yet, we would not like to over-interprete this theory since any SNE in itself 
faces serious interpretation al p roblems, like any other non-linear determin- 
istic Schrodinger equation [18|. Finally, we discussed a novel class of frSNE 
with general complex coupling. Our simulations show that for point-like ob- 
jects the effective (real + induced) gravity is bigger than Newton's by a small 
factor. 



9 



In both our conceptional and numerical analysis we followed restricted 
aims. Without sinking into the interpretational context (wave function col- 
lapse, macroscopic quantum mechanics, quantum gravity, etc.) we wanted to 
point out the hitherto unknown mechanism of gravity induced by the imagi- 
nary mean-field. Our numeric simulations served this aim basically, we have 
not targeted a systematic numeric study of the frSNE. We add, nonetheless, 
the irreversible frSNE is an easy subject of simulation, compared to the irre- 
versible SNE, its numeric study might be an attractive further task in itself 
while, more importantly, it would serve our conceptional understanding a 
possible role of the frSNE in foundations. 
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